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SUMMARY 


The influence of geometry on the isothermal hydrodynamic film separating two rigid 
solids was investigated. The investigation was conducted for a conjunction fully im- 
mersed in lubricant (i. e. , fully flooded). The effect of geometry on the film thickness 
was determined by varying the radius ratio from 1 (a ball -on-plate configuration) to 36 
(a ball in a conforming groove). The dimensionless film thickness was varied from 10”'* 
to 10 . Pressure -viscosity effects were not considered. It was found that the minimum 

film thickness had the same speed, viscosity, and load dependence as Kapitza's classical 
solution. However, the incorporation of the Reynolds boundary conditions resulted in an 
additional geometry effect. That is, the film -thickness equations can be compared as 
follows: 
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where Hq is the dimensionless central (minimum) film thickness, a is the radius 
ratio Ry/R^, e is the film -thickness effect on reduced hydrodynamic lift, and U/W 
is the ratio of dimensionless speed to dimensionless load. With the Reynolds boundary 
conditions the predicted load capacity is 11 to 20 percent greater than if half- 
Sommerfeld boundary conditions .are used. The parabolic approximation results in over- 
estimations of the minimum film thickness of 1. 6 and 0. 7 percent for dimensionless 
minimum film thicknesses of 10 and 10" , respectively. 



INTRODUCTION 


A considerable amount of work (refs. 1 to 6) has been done to develop a minimum- 
film -thickness formula for the classical hydrodynamic point- and line-contact problems. 
Most of the work to date has concentrated on minimum -film -thickness formulations for 
either a ball on a plate or line contact. But the full range of geometries between the two 
extremes has not been adequately studied. Kapitza's film -thickness solution (ref. 1) is 
not limited to a ball -on -plate configuration. However, applying the half -Somm erf eld 
boundary conditions used in Kapitza’s analysis violates flow continuity at the cavitation 
boundary. Consequently, a need exists for a film -thickness formula that is determined 
by more realistic boundary conditions and that applies for a wide range of geometries. 

Work (refs. 7 to 10) has been presented in which the minimum film thickness was 
determined by using an elastohydrodynamic lubrication (EHL) point-contact theory. 

There the influence of contact geometry - as expressed by the ellipticity parameter and 
the effects of the dimensionless speed, load, and material parameters on the minimum 
film thickness - was investigated for both the fully flooded and starved conditions. This 
paper continues the previous work to form a more complete theory. However, to prop- 
erly bring out the physics requires a modification to the analytical approach used in the 
previous work. Here the effect of load capacity has been studied by varying the ratio 
of the transverse radius to the rolling radius (” radius ratio") for two film thicknesses. 
The investigation was conducted for a conjunction fully immersed in lubricant (i. e. , 
fully flooded). The radius ratio was varied from 1 (a ball-on-plate configuration) 
to 36 (a ball in a conforming groove). The dimensionless film thickness was varied from 

e -4 

10 to 10 . Thirteen cases were used in obtaining the film -thickness formula, for 

which pressure and thermal effects were considered to be negligible. Contour and pro- 
file plots of the pressure distribution are also shown. 


SYMBOLS 


a o, a i least -squares coefficients 

D difference, [( H min " H min^ //I *min] x 100, percent 

G material parameter 

H dimensionless film thickness, h/R x 

H q dimensionless central (minimum) film thickness, h Q /R x 

Hq calculated dimensionless central (minimum) film thickness from least -squares 
analysis 

h film thickness, cm 
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hQ central (minimum) film thickness, cm 

L reduced hydrodynamic lift 

m starvation parameter (percent loss in load capacity) 

N direction normal to boundary 

P dimensionless pressure, pR x / i'qU 

O 

p pressure, N/cm 

Q solution to homogeneous Reynolds equation 

R effective radius of curvature, cm 

r radius of curvature, cm 

S separation due to geometry of solids, cm 

U/W ratio of dimensionless speed to dimensionless load 
u average surface velocity in x-direction (u^ +u B )/2, cm/sec 
w load capacity, N 

X dimensionless coordinate, x/R x 

x coordinate along rolling direction, cm 

Y dimensionless coordinate, y/R x 

y coordinate transverse to rolling direction, cm 

a radius ratio, R y /R x 

e film -thickness effect on reduced hydrodynamic lift 
r) dimensionless coordinate, Y/|/2 oHq 

Vq fluid viscosity at standard temperature and pressure, N* sec/cm 
(p Archard-Cowking side-leakage factor, [l + (2/3of)] _1 
X dimensionless coordinate, x/^57 

Subscripts: 

A solid A 

B solid B 

E entrance, or inlet 

x,y coordinate direction 

0 center of contact 

°° infinite domain 



THEORETICAL FORMULATION 


Solution for Central Film Thickness 

The thickness of a hydrodynamic film between two rigid bodies in rolling contact can 
be written as the sum of two terms; that is, 


h = h Q + S(x,y) 


( 1 ) 


where 

hp central (also minimum) film thickness due to hydrodynamic effects 
S(x,y) separation due to geometry of solids 

The separation of two rigid solids (fig. 1(a)) in which the principal axes of inertia of the 
two bodies are parallel can be written as 


S " S Ax + S Bx + S Ay + S By 

where 


S A* = r Ax- ' 

V r L - ** 

S Bx = r Bx - 

V r Bx - ** 

S Ay = r Ay " 

)/ r Ay ' ^ 

s By = r By ' ’ 



A simplifying transformation can be effected by summing the curvatures in the x = 0 
and y = 0 planes. In terms of the effective radius of curvature, 

R x r Ax r Bx 

— = JL+ _L 

R y r Ay r By 
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The resulting equivalent system is shown in figure 1(b). The separation in terms of the 
coordinates and the effective radius of curvature is 


S(x,y) = R x - ^ -x 2 + R y - ^R 2 y - y 2 (2) 

Thus equation (1) is completely determined when the hydrodynamic effects on the central 
film thickness are known. These effects can be determined by applying the conservation 
equations for an incompressible Newtonian fluid under laminar, isothermal, isoviscous, 
steady -state conditions. The following Reynolds equation is obtained: 

JL / h 3 + Jl/h 3 = laivu 2S. (3) 

3x \ ax/ ay \ 3y/ 3x 

where u represents the average surface velocity between the two solids along the rolling 
direction. It is convenient to nondimensionalize with respect to the effective rolling ra- 
dius; that is, 


X s 




(4) 


also 


P = 


v 



where a denotes the "radius ratio. " In terms of these dimensionless variables the 
Reynolds equation becomes 



(5) 


The film -thickness equation in dimensionless form is 



2 2 

For situations in which X « 1, and (Y /a) « 1, it is convenient to expand H in a 
two-dimensional Taylor series to give 
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( 7 ) 


■y 2 y2 
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This is called the parabolic -approximation. The analysis to follow uses both forms of 
H, and comparisons are made. 

The solution of the Reynolds equation (eq. (5)) is known to consist of a homogeneous 
solution and a particular solution; that is, 


P = Pp + P h (8) 

for which is a solution to the homogeneous equation and satisfies the condition that 
P^ = -Pp at the boundaries: 



( 9 ) 


For the parabolic approximation, the particular solution for the pressure is simply pro- 
portional to X/H^; that is, 


where 


P - 


IT 



In the preceding equation <p is the side -leakage factor established by Ar chard and 
Cowking (ref. 11) and can be verified by inserting P p back into equation (5). If we de- 
fine P^(X, Y) = 4 <^Q(X, Y), by using equation (8), we can express the full solution as 

P-4?(:£h-q\ (10) 

Vh 2 1 

In general, the homogeneous solution P^ is an unknown function of X and Y. Conse- 
quently, the pressure distribution must be determined numerically. The pressure dis- 
tribution as given in equation (10) can be used, however, in relating the hydrodynamic 
effects (i. e. , load, speed, and viscosity) to the central film thickness. First, the load 
capacity and the pressure distribution are related through the following equation: 
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w = ff p dx dy 


(ID 


or in dimensionless form (using eq. (4)) 

w = i/ q uR x //P dX dY 

Substituting equation (10) into this expression gives 



For the parabolic film assumption, the central film thickness can be isolated from the 
integrand by defining the following transformation: 


X = |/2Hq x 
Y = >/2aH 0 rj 


If we assume the homogeneous solution to transform in the same manner as the par- 
ticular solution, we obtain 


w = 8^ 0 uR x */_ 



* + Q(x, V) 

(l + x 2 + » 2 ) 2 


d\ dri 


( 12 ) 


Kapitza refers to this integral as the reduced hydrodynamic lift L. Thus, 

d X dv (13) 

The reduced hydrodynamic lift in Kapitza's analysis was determined to be w/2 by as- 
suming Q = 0 and integrating over the half -space of positive pressures. For Reynolds 
boundary conditions, the limits depend on the shape of the cavitation boundary and hence 
the geometry. Consequently, we seek an additional geometry effect in Kapitza's solution 
Equation (12) enables us to determine the central (minimum) film thickness as a function 
of the load, speed, geometry, and fluid viscosity; that is, 


L s 


— =*- 
JJ Ld + x 2 + n 


2) 2 


+ Q(x,v) 
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For the parabolic approximation, we need to determine L only as a function of the 
geometry; that is, 

(x 2 « (1/2H 0 ) 

L = L(a) if < 

|/ « (1/2H 0 ) 

This will be determined numerically. If, on the other hand, the film thickness is large 
enough that these inequalities cannot be satisfied throughout most of the domain, the 
exact film -thickness equation (eq. (6)) must be used. The integrand of equation (13) thus 
becomes a function of the central film thickness. Consequently, L = L(a,Hg), resulting 
in a transcendental equation for Hq. 


Boundary Conditions 

Earlier theories (ref. 12) assumed the pressure to be ambient or zero at the point 
of closest approach. This resulted in an antisymmetric solution with respect to X 
(fig. 2(a)). In actuality, the lubricant is unable to sustain the negative pressures pre- 
dicted by the full solution. A simple approach taken by Kapitza (ref. 1) was to ignore 
the negative pressures, that is, to employ the half -Somm erf eld boundary conditions. 

This solution (fig. 2(b)) has been used to get a reasonable estimate for the load capacity. 
However, Kapitza's solution does not satisfy continuity conditions at the exit (cavitation) 
boundary; that is, the pressure gradient normal to the cavitation boundary must be zero. 
To insist on P = (9P/9N) = 0 at X = 0 would be overspecifying the problem mathe- 
matically. However, we can insist on P = (9P/9N) = 0 at the cavitation boundary 
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(i. e. , Reynolds boundary conditions). The general solution will then appear as in fig- 
ure 2(c). For film thicknesses of about 10 to 10 centimeter, Dowson (ref. 13) has 
shown that other boundary requirements are needed. For that investigation, the cav- 
itation boundary as determined by the Reynolds boundary conditions did not coincide with 
the cavitation boundary observed experimentally. Taylor (ref. 14) has summarized sev- 
eral of the boundary requirements that have evolved as a result of Dowson's work. 

Among the boundary requirements discussed, the Reynolds boundary conditions were 
chosen as the most appropriate for the conditions of load and speed in this investigation. 


NUMERICAL ANALYSIS 

A pressure distribution satisfying the Reynolds equation (eq. (3)) was determined 
numerically for a given speed, viscosity, geometry, and film thickness. The numerical 
solution was achieved by using the Gauss -Seidel iterative method with overrelaxation. 
The parameter $ = (ref. 15) was introduced to help the relaxation process. A 

computer program is given in the appendix. 


Nodal Structure 

A variable -mesh nodal structure (fig. 3) was used to provide close spacing in and 
around the pressure peak. This helped to minimize the errors that can occur because of 
large gradients in the high-pressure region. The grid spacing used in terms of the co-‘ 
ordinates X and Y was varied depending on anticipated pressure distribution. That is, 
for a very highly peaked and localized pressure distribution, the fine mesh spacing was 
about 0. 002 and the coarse mesh spacing was about 0. 1. For a relatively flat pressure 
distribution, the fine mesh spacing was about 0. 005 and the coarse mesh spacing was 
about 0. 13. 


Integration Domain 

The size of the conjunction, or the integration domain, is determined so as to make 
the contact fully flooded or as close to that as practical. From Dalmaz (ref. 3), a fully 
flooded condition for a ball -on -plate configuration would have an inlet domain defined by 
Xg = -1 and Yg = ±1. By using the pressures from Kapitza's classical theory, the loss 
of load capacity resulting from using a finite inlet rather than the semi -infinite inlet used 
by Kapitza can be estimated. For this purpose let m represent the percent loss in load 
capacity; that is, 


9 



m = 



(17) 


Thus as m approaches 100, the inlet becomes severely starved. If m approaches 
zero, the inlet is considered fully flooded. For the ball -on-plate inlet domain, m is 
calculated to be 1. 61 percent. Since this represents a negligible loss in load capacity, 
we chose to retain the concept that a fully flooded condition exists in this case. Hence- 
forth, if m < 1. 61 percent, the inlet is considered as fully flooded. According to this 
criterion, the inlets for all the geometries considered in this investigation were fully 
flooded (table I). The exit boundary was determined so as to allow for a fully developed 
cavitation boundary. 

In our effort to achieve a fully flooded condition, we recognize the fact that the 
Reynolds equation loses some of its validity at large distances from the point of minimum 
film thickness. Dowson (ref. 16) has pointed out that the errors involved in using this 
equation to determine the buildup of pressure in such regions are negligible: The pre- 
dicted pressures are themselves so very much smaller than the effective load -carrying 
pressures in the region of closest approach of the solids. 


Film Thickness (Parabolic Film Approximation) 

Once the integration domain has been established, the film -thickness equation can 
be determined numerically. First, the load capacity w is obtained from the numer- 
ically determined pressure distribution. Then inserting the value of w into equation (14) 
allows us to solve for L for various geometries. 

From the data in table II, a curve fit can be effected as a function of geometry. 
Studying figure 4 and trying several appropriate functional forms to curve -fit the data 
showed that the following equation represented the data best: 

L = a x tan" 1 - + a Q (18) 

2 


The values of aQ and a^ were determined to be 


a Q = 1.683^ 

► 

a. = 0. 131 

J 


(19) 


o 9 

The coefficient of determination r (ref. 17) was 0. 96 for this fit. The value of r re- 
flects the fit of the data to the resulting equation: 1 being a perfect fit, and 0 being the 
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worst possible fit. Inserting equation (18) into equation (16) gives for the calculated min- 
inum film thickness 


H min 


= 128a 


22 (0. 131 tan' 1 2 + 1.683 
W \ 2 


( 20 ) 


For comparison of the calculated film thickness with the actual input film thickness, it is 
convenient to define the differences (in percent) as 


D = 


( H min " H min^ 


x 100 


H 


min 


The difference ranged from -2. 14 to 1. 35 percent (table n). 


( 21 ) 


Film Thickness (Exact Equation) 

The results of using the exact film -thickness equation (eq. (6)) rather than the para- 
bolic approximation (eq. (7)) can be compared from the values of L given in table n. 

The values of L determined through the exact film -thickness analysis are always re- 
duced by a constant factor. Thus the dependence on geometry and film thickness for the 
reduced hydrodynamic lift can be separated as follows: 

L(a,H 0 ) =<r(H 0 )L p (a) (22) 

where Lp is determined by using the parabolic approximation. The values of £(Hq) are 

given in table in. Using the parabolic approximation gives errors of 0. 7 and 1. 6 percent 

-5 -4 

for film thicknesses of 10 and 10 , respectively. Although these errors are quite 

small, they would probably be larger for thicker films: As the film thickness is in- 
creased, the pressure distribution spreads out more evenly (e. g. , fig. 5). Thus the 
pressures far from the point of contact, where the parabolic assumption is no longer 
valid, contribute more to the load capacity than if the pressure distribution had been very 
localized. 

The minimum film thickness can now be written as 


H. 

mm 


= 128a 


[ (o. 131 tan' 1 -+ 1. 683 
W \ 2 


(23) 


The difference ranges from 2. 00 to -2. 11 percent for the exact film -thickness results in 
table n. 
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DISCUSSION OF RESULTS 


Comparison of Theories 

The minimum -film -thickness equation derived by Kapitza using half -Somm erf eld 
boundary conditions and assuming the parabolic approximation is 


H. = 128 a m - 

mm l W2 


By equating equation (16) with this expression for a given speed parameter, the load ca- 
pacity for the two theories can be compared; that is, 



(25) 


The effect of the geometry on the reduced hydrodynamic lift is shown in figure 4. The 
figure shows that L, and hence the load capacity (eq. (25)), is 11 to 20 percent greater 
than that predicted by Kapitza. The least difference occurs for a ball -on -plate contact. 
As a is increased, the difference in load capacities approaches a constant 20 percent. 
The alteration of the pressure distribution due to the Reynolds boundary conditions at the 
cavitation boundary is responsible for this geometry effect. Figure 6 is a three- 
dimensional representation of a pressure distribution for a of 1. 00 and 36. 54 and il- 
lustrates the shape of the cavitation boundary. As a becomes large, the cavitation 
boundary tends to straighten out, accompanied by decreasing changes in L. The scale 
along the y-axis in figure 6(a) has been magnified about 3 times to improve the resolu- 
tion. Consequently the differences in the shapes of the cavitation boundary are actually 
subdued as they are presented. 

The two analyses resulted in the same exponent of 2 for U/W. Dalmaz(ref. 3), using 
the Reynolds boundary conditions for a ball-on-plate configuration, reported an exponent 
of 1. 77. The lower exponent appears to be due to starvation effects resulting from the 
inlet condition in both the analytical and experimental results. This is illustrated more 
clearly by comparing the results of applying the inlet condition of Dalmaz and the inlet 
condition in this study, using Kapitza's classical theory for both. A fixed oil film thick- 
ness at the inlet that is independent of the minimum film thickness is used in Dalmaz’s 
analysis and in this study as well. But the inlet oil film thickness for the analysis of ref- 
erence 3 and experimental work (given in ref. 18) was roughly 5 percent of that used 
here. The effect that this has on the exponent of U/W can be seen by calculating the 
load capacity for two different film thicknesses while keeping the other hydrodynamic 
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variables constant. Dimensionless minimum -film -thickness values of 10 and 10 re- 
sult in load capacities of 0. 046 and 0. 151 newton, respectively. This would yield an ex- 

-3 

ponent of 1. 94. For thicker films (i. e. , Hq = 10 ) the starvation effects become even 

more pronounced and the U/W exponent is driven down in value to 1. 84. For com- 
parison, from the data in table I for fully flooded conditions, we calculate the exponent 
to be 1. 98. 

Figure 7 compares experimental data (ref. 3) taken under lightly loaded (rigid con- 
tacts), isoviscous conditions for pure sliding of a ball on a plate with the corresponding 
theoretical results. The parameter grouping in the figure is that used by Dalmaz (ref. 3) 
and first introduced by Thorp and Gohar (ref. 4). The theoretical results of this paper 
are in excellent agreement with the data for the lower half range of Hq/WG. For the 
reasons previously explained, the agreement of these results with the experimental re- 
sults of reference 3 begins to diverge for the upper range of Hq/WG. Further, using 
comparable inlet conditions for both theory and experiment provides better agreement 
between the two for the upper range of Hq/WG, as the theoretical line of Dalmaz attests 
in figure 7. For comparison, the theory by Kapitza has been included. 


Contour Plots 

( Isobar plots for three radius ratios (i. e. , a of 25. 29, 8. 30, and 1. 00) are shown in 
figure 8. The contours were generated by means of a contour -plotting subroutine and 
displayed with a Cal comp plotter. The contours belong to the family of curves defined by 
equation (10). The center of contact is represented by the asterisk. The pressure peak 
builds up in the entrance region, which is located to the left of the center of contact and 
is indicated by the +. Since the isobars in each case are evenly spaced, the pressure 
gradients can be easily depicted. Note that, as the radius ratio increases, the steeper 
pressure gradients are predominantly along the rolling direction. This implies that the 
amount of side leakage decreases as a increases. A decrease in side leakage is re- 
flected in an increase in the value of cp. For line contact cp = 1 and for the largest 
value of a in this investigation (a = 36. 54) cp = 0. 998. 


Pressure Profiles 

Pressure profiles across the center of the conjunction and in the direction of roUing 
are shown in figure 9 for three radius ratios (i. e. , or of 36. 54, 2. 84, and 1. 00). The 
pressures were generated for a constant dimensionless film thickness (Hq = 10 ~^). The 
locations of the pressure peaks were not altered by the fact that Reynolds boundary con- 
ditions were used rather than half -Somm erf eld boundary conditions. The locations of the 
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pressure peaks from reference 1 are determined by setting (3P/3X) = (SP/3Y) = 0 and 
solving for X and Y as follows: 





(26) 


Equation (26) shows that the geometry does not affect the pressure peak location, as can 
be verified by the numerically determined curves of figure 9. However, the magnitude 
of the pressure peaks will, according to equation (10), be modulated by the Archard- 
Cowking side-leakage factor <p, which is a function of the geometry. The influence of 
the side -leakage factor on the pressure distribution for several geometries is shown in 
figure 9. 


CONCLUDING REMARKS 


The influence of geometry on the isothermal hydrodynamic film separating two rigid 
solids was investigated. The investigation was conducted for a conjunction fully im- 
mersed in lubricant (i. e. , fully flooded). The effect of the geometry on the film thick- 
ness was determined by varying the radius ratio from 1 (a ball -on -plate configuration) to 
36 (a ball in a conforming groove). The dimensionless film thickness was varied from 
Hr 5 to lO -4 . Pressure -viscosity effects were not considered. It was found that the min- 
inum film thickness had the same speed, viscosity, and load dependence as Kapitza's 
classical solution. However, the incorporation of the Reynolds boundary conditions re- 
sulted in an additional geometry effect. That is, the film -thickness equations can be 
compared as follows: 


H Q = 128a 


e(H n ) SdH (o. 131 tan -1 - + 1. 683 
u W \ 2 



{ Full circular film; 

Reynolds boundary conditions 


Hq = 128a 


^ ( 0 . 131 tan" 1 -+ 1.683 
W \ 2 


Parabolic approximation; 
Reynolds boundary conditions 



Parabolic approximation; 
half -Sommerfeld boundary 
conditions 
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where Hq is the dimensionless central (minimum) film thickness; a is the radius ratio 

R /R ; e is the film -thickness effect on reduced hydrodynamic lift; cp is the Archard- 
y x 

Cowking side-leakage factor; and U/W is the ratio of dimensionless speed to dimension- 
less load. The Reynolds boundary conditions resulted in the predicted load capacity 
being 11 to 20 percent greater than if half -Sommerf eld boundary conditions were used. 

The parabolic approximation resulted in overestimations of the minimum film thick- 
ness of about 1. 6 and 0.7 percent for calculated dimensionless film thicknesses of 10"^ 

_5 

and 10 , respectively. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 12, 1978, 

505-04. 
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APPENDIX - COMPUTER PROGRAM 


IMPLICIT REAL*8 (A-H, 0-Z) 

REAL*4 XM IN, XLSN,YMIN , XE (60) 

RE AL*4 XG (56) , YG ( 56) , YGC ( 1 1 9) , DtJMM (56 , 1 1 9 , 2) 

RE AL*4 FN,U5,WS,3S,XKS,XMAX,XCENTS,TCENTS,YMAX,PRMXS,XPK 
REAL*4 XN0DFS,XD3L,YDEL,HMINS,PMIN,HMAX 
NAMELIST/ PADIUS/RAX, RAY, PBX,RBY,RX,RY,U A, UB, PHEE 
NAMELIST/OUTPUT/P, PBAR ,PCHECK,SUM3 ,H0 ,HOBAR, OOW 
DIMENSION P2 (56,50,2) 

DIMENSION DENS (336 0) , VIS (3 360), XMU (33 60) ,ZPP (3360) , PHI (336 0 1 
DIMENSION A( 3360) , B( 3360) ,C(3 360) , DLZ (3 3 6 0) , XL (33 6 0) , XM (3 3 5 0 ) 
DIMENSION S(3 36 0) , H (3 360) ,PPSV (3360) , PR (3360) 

DIMENSION IN (3) 

COMMON/CTRINF/FN,U3,W5,GS, XM A X , XCENT S, Y CENTS , Y M AX , PR « XS , X PK 
COMMON /C TR/XNODFS, XD E L , Y DEL, H MI NS , PMI N, HM AX 
COMMON/CTRINT/NCGF 

INPUT 

MPARA= 1 
N V ISC=0 
MWT= 1 
C 

C MPARA= 1 : Tii E PAR A30LIC APPROXIMATION MILL BE USED 

C NVISC=0: THE ISO-VISCOUS SOLUTION IS DETERMINED 

C MWT= 1 : PARTIALLY CONVERGED SOLUTION IS STORED ON TAPE 

TN CASE OF CRASH 
C 

ROSM=1.0D0 
IDE= 1 
M A U R = 0 
JE NN=0 
ZZZ= 1 . 0D0 

PI=3. UM 59255359D0 

PIVAS=4358.7D0 

VISO=4. 1 1D-6 

VA=0.0D0 

VB = 0 . 0 DO 

Z=.67D0 

ELPHA=5. 82744D-6 
EET A =1,683430-5 
VISE=. 000000006.31 DO 
OR F= 1 , 9D0 
NX=56 
NY =60 

ZD1=9. 799652 DO 
ZE 1=490. 1 96078D0 
ZE2=ZD1 
ZC1=10. 000D0 
ZC2=500- 0D0 
JZON= 1 0 
J MI N= 1 0 
NCGF= 1 0 
NZON2=44 
XCENT= 1 . OODO 
YCENT= 1 . OODO 
XNODF= . 9134000 
M= 1 


16 



c 

C M=-1: READS IN FORMATTED INPUT AND BYPASSES THE 

INITIAL CALCULATION OF 
C THE PRESSURE DIS T R TOUT ION 

C 

C M = 0: READS IN H EX A DECIMAL INPUT AND BYPASSES THE 

INITIAL CALCULATION OF 
C THE PRESSURE DISTRIBUTION 

C 

C M= + 1: CALCULATES AN INITIAL GUESS 

C 

R A X= 1 . 11125D0 
RAY=RAX 
RBX= 1 . D12 
RBY = - 1.D12 
UA= 1 0. 0D0 

975 CONTINUE 
DEL1=. IDO 
MAXS2=1 
D2L3=. IDO 
HM IN = 1 0 • 0 DO 
HMAX=10.0 

PR S VM X = . ID-18 
PRMX=- ID-18 

976 CONTINUE 
UB=UA 

NX Y = NX*N Y 
NX 1 = NX- 1 
NY 1=2*NY 
NYM1=NY1 -1 
IN Y = 2* NX 
7.C 1 S = ZC 1 
ZC2S = ZC 2 
ZE1S=ZE1 
ZE2S=ZE2 
FN=FLOAT (IDE) 

IN ( 1) =IDE 
IN (2) = NY 
IN (3) =NCGF 
C 

C CURVATURE SUM AND DIFFERENCE 

C 

RX= (1./RAX) (1./RBX) 

RY= (1. /RAY) + ( 1 . /R 3 Y) 

RHO=RX+RY 
G A M M A = (RX-RY) /RHO 
C 

XCENTS=XCENT/RX 

YCENTS=YCENT/RX 

XNODFS=XNODF/PX 

YM=YCENT 

XN=XCENT 

C 

C DIMENSIONLESS PARAMETER GROUPING 

C 

SMU = . 5* (OA+UB) 

SMV = . 5* (VA + VB) 

V=SQRT (SMU**2+SMV**2) 

THETA= ATAN (3MV/SMU) 

EP = VISO* V*RX 
U= VISG*RX*V/EP 
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us=u 

H0=1.D-4 , 

PHEE = 1 • / ( 1 .♦ (2.*RY)/(3.*RX)) 

Q1 = VISE/VISO 
Q2=EP/19608. 5268 
Q5=ELPHA*EP 
Q6=BETA*EP 

4 CONTINUE 
ZC=ZC1 
ZCS=ZC 
SA=1./ZC 
SB= 1.0 /ZD 1 

INITIAL GUESS 

IF (M) 5,6,7 

5 CONTINUE 

READ (5,7 00) ( ( (P2 (I,J,K) , 1= 1, NX) , J= 1 , NY ) ,K=1,2) 

700 FORMAT (8D10. 5) 

GO TO 7 

6 CONTINUE 

READ (5) ( ( P 2 ( I , J , 1) ,1 = 1, NX) ,J = 1,NY) 

READ (5) ( (P 2 (I, J, 2) , 1=1, NX) , J=1,NY) 

REWIND 8 . 

7 CONTINUE 
XMIN=XCENTS-. 20 
Y=0. DO 

DO 13 J= 1 , NY 
X=0. DO 

DO 12 1=1, NX 
ZD= ZD 1 

N=I+ (J-1) *NX 
JTWIN=NY1-J 
IF ( J . GE. JZON) ZC=ZC2 
YGC ( J) = Y/RX 

YGC (JTWIN) =2. +YCENTS- YGC ( J) 

YMIN= YCENTS-. 16 
YHAX=YCENTS+. 16 
XG (I) = X/RX 

IF (YGC (JTWIN) .GE.YMAX) JMAX=JTWIN-1 
IF (JTWIN. EQ. JMAX .AND. XG ( I) . L E. XCSN TS ) LU = N 

IP (XG (I) . LE. XKIN .AND. J. EQ. NY) K 1=N 

IF(XMIN.EQ.O.) K1=NX* (NY-1) +2 

IF (J. EQ. J MIN .AND. X G ( I) . LE. XCENTS) LM AX = N 

NVMF=NCGF* (J-1) *NX 

NVMB=NVMF«-NZON2 

IF (N. GE. NVMP .AND. N. LT. NVMB) ZD=ZE1 

IF(N.GE.NVMB) ZD=ZS2 

SQX= ( (X-XN) /RX) **2 

SQY= ( (Y-YM)/RX) **2 

IF (MPARA .EQ. 1) GO TO 75 

S (N) = (1. /RX) ♦ (1. /R Y) -DSQRT( (1./RX) **2-S QX) - 0SQRT( (1./RY) **2-SQY) 
IF (S(N) .GT. 1./RX) S (N) =1 ./RX 
GO TO 8 

75 S(N) = (RX*SQX«-RY*SQY)/2. 

IF (S ( N) .GT. 0.5/RX) S (N) =0. 5/RX 

8 CONTINUE 

IF (M-EQ. 1) 30 TO 9 
PR(N) =P2(I,J, 1) 

GO TO 10 



non ono oon 


9 CONTINUE 
PR (N) = - ID-8 

10 CONTINUE 

IF (J.EQ.1 .DR. I. EQ. 1 .OR. I.EQ.NX) PR(N)=0.0D0 
IF (PR (N) . GT. 0 . ) 3D TD 11 
PR (N) =0. DO 

11 PRSV(N) =PF (N) 

IF(D.EQ.NY) YG(I)=PRtN) 

IF (PR (N) .LT. PRSVMX) GO TO 12 
PRS VMX= PR (N) 

NHOLD=N 

12 X = X+ 1 . /ZD 

13 Y=Y+1./ZC 
XMAX=XG (NX) 

YDEL=YMAX-YMIN 
XDEL= YDEL 

WRITE (6, 900) (PR(N) ,N=1,NXY) 

900 FORMAT (1H ,2HPR/70(1H , 10 D 1 3. 5/) ) 

WRITE (6, 1000) (S(N) ,N=1,NXY) 

1000 FORMAT ( 1H ,1HS/7D(1H ,10013.5/)) 

FILM THICKNESS 

14 CONTINUE 

DO 19 J= 1 , NY 
DO 18 1=1, NX 
N3 = I+ (J- 1) *N X 
H (N3) =HO+RX*S (N3) 

IF (H ( N 3) .GT. HMIN) 3D TO 17 
H MI N= H (N3) 

HMINS=HMIN 
NS A VE= N3 

17 CONTINUE 

PHI (N3) = PR (N3) * (H (N3) *+1.5) 

18 CONTINUE 

19 CONTINUE 

WRITE (6,1 100) NS A V E, HMIN 
1100 FORM AT ( 1 H , 6 H N S A V E = , I 1 0, 1 OX, 5 HHMI N= , D 1 6 . 5 ) 

WRITE (6,1200) PRS VKX , NHOLD 
1200 FORMAT(1H ,8HPRSV'1X=,D16.5, 1 0 X, 7H NHOLD= , 1 1 0) 

20 CONTINUE 
WRITE(6, RADIUS) 

INITIAL VISCOSITY AND DENSITY CALCULATION 
DO 21 N = 1 , NX Y 

DENS (N) = 1. *■ (Q5*P8 (N) ) / (I. +Q6*PR (N) ) 

VIS(N) =Q1**(1.- (1.+Q2*PR(N) ) **Z) 

XMU (N) =DENS (N) /VIS (N) 

21 CONTINUE 

RELAXATION CDEFFIENTS A , D , C, D , L, A ND M 

22 CONTINUE 
SUM^O. 0D0 
ZC=ZC1 
ZF=?.C1 

DO 24 J=2 , NY 
DO 23 1=2, NX1 
ZD=ZD1 
ZE=ZD 
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N = I + (J-1) *NX 
NVMF=NCGF+ (J-1) *NX 
NVMB=NVMF+NZ3N2 
IF (J. RE. JZON) ZF=ZC2 
IF (J.GT. JZON) ZC=ZC2 

IF (N. GE. NVHP .AND. N. LT. NVMB) ZE=ZE1 
IF (N. GT. N VMF .AND. N. LE. NVMB) ZD=ZE1 
IF (N. GE. NVMB) ZE=ZE2 
IF (N.GT. NVM3) ZD = ZE2 
V 1 = ( (ZE*ZD) / (ZE+ZD) ) **2. 

V2= ( (ZE) ** 2 . ) * (ZD-ZE) /(ZD+ZE) 

V3= (ZE**3. ) * (ZE+2. *ZD) /( (ZE+ZD) ) **2. 

V4= (ZD**3. ) * (2. ♦ZE+ZD) /( (ZE+ZD) ) **2. 

V5= (ZD/ZE) **2.*V2 
C2D=ZC**2-ZF**2 
C3P1=ZC+2.*ZF 
C3P2=2. *ZC+ZF 
CF= (ZF/ (ZC+ZF ) ) **2 
CF 1 = (1. / (ZC+ZF ) ) **2 
CP= (ZC + ZF ) **2 
Z1=. 5 * 7. D 
Z2=ZC/ (ZC+ZF) 

Z3=Z1**2 
ZU=Z2**2 
N 1=N+ 1 
N2=N-1 
N3=N+NX 
N4=N-NX 

IF (J.EQ.NY) N3=N4 
YO=XMU (N) 

Y 1 =XM U (Nl) 

Y2 = XMIJ (N2) 

Y3=XMU ( N 3 ) 

Y4=XMU (N4 ) 

Y5=H (N 1) 

Y6 = H (N2) 

Y7=H (N) 

Y8 = H (N3) 

Y9=H ( N 4 ) 

Y10=Y1*DSQRT (Y5) 

Y 1 1=Y2*DSQRT (Y6) 

Y12=Y3*DSQRT (Y8) 

Y 13=Y4*DSQRT (Y9) 

Y14=YO*DSQRT (Y7) 

A (N) =V1*Y2+V2*YD+V3*Y1 
B (N) =Z4*(Y3*ZF**2-C2D*YO+ZC*C3P1*Y4) 

C (N) =V4*Y2-V5*Y0+ V 1*Y 1 

DLZ (N) =CF* (ZP*C3P2*Y3+C2D*YO+ZC**2*Y4) 

XL 1 = Y2*4 . * 7 . 3 - (ZD-ZE) **2*Y0+Y1*ZE**2+CP1* (Y3*ZF**2*CP-:2D**2* Y0+ZC**2*CP*Y4) 
XL2=1.5/(Y7**1.5) 

XL 3 = Y 1 1*(V4*Y6-Y7*4.*Z3 + Y5*V1)+Y14*(-V5*Y6+ (ZD-ZE) **2*Y7+Y5* V2) + 

1 Y10* (V 1* Y6-ZE**2*Y7+V3*Y5) 

XL4=CF1* (Y12*ZF**2*(ZF*C3P2*Y8-CP*Y7+Y9*ZC**2) ♦ 

1 Y13*ZC**2* (Y3*ZF**2-CP*Y7+ZC*C3P1*Y9) + 

1 Y14*C2D* (ZF**2*Y8+C2D*Y7-ZC**2*Y9) ) 

XL(N)=XL1+XL2*(XL3+XL4) 

XH1 = 12.*0/(Y7**1.5) 

XH2=( (ZE**2. *DENS (Nl) *Y5/ (ZD + ZE) + (ZD-ZE) * DENS (N) *Y7- 
1 ZD**2. *DE NS (N2 3) *Y6/ (ZD + ZE) ) *DCOS (THETA) ) 

XH 3- (Y8*ZF**2*DENS (N3) +C2D*DENS (N) * Y7-Y9 *ZC* *2*DENS (N4) )* DSTN (THETA) 
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XH (N) =XH1* (XH2 + Z2/ZC*X«3) 

23 CONTINUE 

24 CONTINUE 

25 SUM=0.0D0 

RELAXATION FORMULA 

DO 30 J = 2 ,N Y 
MEND=0 

DO 29 1 = 2, NX1 

MN = I + (J-1) *N X 

MNA=NN+1 

MNB = MN- 1 

SNC=MN+NX 

HND=MN-NX 

IF(J.EQ. NY) MNC=MND 

ZPR (MN) =PHI ( !1 N) -0RF* (PHI (MN)M*M(MN)-A(MN) *PHI (MN A) -B ( MN) * PH I (H N D) - 
1 -C(MN) *PHI (MNB) -DLZ (MN) *PHI (NNC) ) /XL(MN) ) 

26 CONTINUE 

IF (ZPR (MN) .LE.O.) GO TO 27 
Y18 = (ZPR (PIN) - PHI (Mil) ) /ZPR (MN) 

SUM=SUM+DABS ( Y18) 

GO TO 28 

27 ZPR (MN) =0. ID-20 
MEND=HEND*-1 

IF (MEND .E3. 1) XE (J) =XG (I) 

28 PHI (MN) =ZPR(MN) 

29 CONTINUE 

30 CONTINUE 

31 CONTINUE 
MAUR=MAUR+ 1 

41 CONTINUE 

IF (MWT.EQ.O) GO TO 66 

65 CONTINUE 

WRITE (8) ( ( P 2 (I,J,1) ,1 = 1, NX) , J = 1 , N Y) 

WRITE (8) ((P2(I,J,2) ,1 = 1, NX) ,J=1,NY) V 

REWIND 8 

66 CONTINUE 

IF (SUM.LT.DEL3) GO TO 311 
GO TO 25 

VISCOSITY AND DENSITY ITERATION 

311 SDM2=0. 0D0 

WRITE (6,480) MAUR.SUM 
480 FORMAT (1H , 5 HM A UR =, 1 8 , 5 X, 6H SUM= , 18) 

MAUR=0 
PRMX=. ID-13 
DO 315 J=2 ,NY 
DO 316 1=2, NX1 
N=I+ (J-1) * NX 

PR (N) = (PHI (N) /(R(N) *«■ 1.5) ♦ (ROSM-1.) *PR(N) ) /ROSM 
IF (PR (N) . LT. PRHX) GO TO 33 
PRMX = PR (N) 

PRMXS=PRMX 
NHOLD=N 
33 CONTINUE 

IF(NVISC.EQ.O) 30 TO 316 
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DE NSN= 1 . ♦ (Q5*PR(N> )/(1.*Q6*PH(5J)) 
VISN=Ql** (1.- (1. + 22*PR (N) ) **Z) ’ 
XMUN=DENSN/VISN 
Y99= (XBUN-XMU (N) ) /XMUN 
SUM2=SUM2+DABS (T99) 

DENS (N) =D£NSN 
VIS (N) =VISN 
316 XMU (N) = XMUN 
315 CONTINUE 
MN 1=1 

459 CONTINUE 

JENN=JENN+ 1 
WRITE (6,810) SUM2 
810 FORMAT (1(1 , 3 HSUM2 1 , D1 G. 5) 

IF (SUM2.LT. DELI) ' 30 TO 32 
IF (MAXS2. GT. 400) GO TO 63 
MAXS2=MAXS2* 1 
GO TO 22 

APPLIED NORMAL LOAD 

32 CONTINUE 
MAXS2= 1 

WRITE (6,481) JF.NN 

481 FORMAT ( 1H , 9H JENN IFER =, 18 ) 

JE NN=0 

M A UR = 0 
QU4=0. ODO 
PSUM=0. ODO 
IZON=1 
ZD=ZD1 
PBAR=0. ODO 
PCHECK=0. ODO 
MS=2 
MF=NCGF 
GO TO 483 

482 ZD=ZE1 
I ZON= 2 
MS = MF 

MF=MS+ NZON2 
QU4=0.0D0 
PSUB=0. 0 DO 
GO TO 483 

485 ZD=ZE2 
IZON=3 
MS = MF 
MF=NX 
QU4 = 0. 0 DO 
PSUM=0. ODO 

483 CONTINUE 

DO 38 1 = M S , M F 
I M= MOD (1,2) 

IMS=MOD (MS, 2) 

IF (IMS. EQ. 1) GO TO 486 
QU 1= 1 . ODO 

IF(IM.EQ.I) QU1=2. ODO 
GO TO 487 

486 QU 1=2. ODO 

IF (IM . EQ. 1 ) 2U1=1 . ODO 

487 CONTINUE 
QU5=1 . ODO 



QU3=0. ODO 

IF(I.EQ.NS .OR. I.EQ.MF) QU1=.5D0 
IF (I. EQ. MS .OR. I.EQ. OF) QU5=.5D0 
JINTG=1 
ZC=ZC1 

NS = 2 
NF= JZO N 
GO TO 35 

34 JI NTG=2 
ZC=ZC2 
NS=NF 
NF = NY 

35 DO 37 J= NS , N F 
Hal* (J-1) *NX 
JN=MOD (J,2) 

JNS=MOD (NS, 2) 

IF (JNS. EQ. 1) GO ID 40 
IF(JN.EQ.O) Q V 1 =2 . ODO 
IF (JN. EQ. 1 ) QV1=4.0D0 
IF ( J . EQ. N S .OR. O.EQ.NF) QVl=1.OD0 
GO TO 37 

40 IF (JN . EQ. 0) QV1 =4 . ODO 
IF ( JN. EQ. 1 ) QV 1 =2. ODO 
IF ( J. EQ . N S .OR. J. EQ. NF) QV1 = 1.0D0 

37 QU3=QU3+PR(N)*QV1/ZC 

IF (JINTG.EQ.1) GO TO 34 
PSUB= PSUB+QU3*QU5/3. 0 

38 QU4=QU4*QU3*QU1/3.0 
PBAR=4.*EP/(RX**2) *QU4/ (3 . *ZD) *PBAR 
PCHECK = 2. *EP/ (RX* * 2) * PSUM/ZD+PCHECK 
IF (IZON. EQ.2) GO TO 485 

IF (IZON. EQ. 3) GO TO 36 / 

GO TO 482 

36 CONTINUE 
UOW=U*VISO/(PBAR*RX) 

WRITE (6, OUTPUT) 

WS=PBAR 

39 CONTINUE 

IF (MPARA.EQ.O) GO TO 43 
PRMX=. ID-13 
PRSVMX=PRMX 
H M AX= 1 0. 0 DO 
HMIN= 1 0 . 0 DO 
MPARA=0 
GO TO 976 
43 CONTINUE 

WRITE (6, 1600) NS A V E, H HI N 
1600 FO RBAT (1 H , 6 H NS AV E= , 1 1 0 , 1 0X , 5 HHMIN=, D 1 6. 5) 

WRITE(6, 1700) (PR (N) ,N=1,NXY) 

1700 FORMAT (1 H ,2 flPK/7 0 ( 1 H ,10D13.5/)) 

WRITE (6, 1800) (PHI (N) ,N=1 ,NXY) 

1800 FORMAT ( 1 H ,3HPHI/70(1H ,10013.5/)) 

WRITE (6, 1900) (H(N) , N = 1 ,NXY) 

WRITE (6,2000) (S ( N) , N= 1 , NX Y) 

1900 FORMAT ( 1 H ,1HH/70(1H ,10013.5/)) 

2000 FORMAT (1H ,1HS/70(1H ,10D13.5/)) 

WRITE (6, 2100) (DENS (N) ,N=1,NXY) 

2100 FORMAT (1 H , 7 HDEN3 IT Y/70 ( 1 H ,10D13. 5/)) 

WRITE (5 ,2200) (VIS (N) , N=1 ,NXY) 

2200 FORMAT ( 1H ,9HVISCOSITY/70 ( 1H ,10D13.5/)) 

WRITE (6,2300) (A ( N) , N= 1 , NXY) 
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2300 

2400 

2500 

2600 

2700 

2800 

C 

53 


57 


58 


54 

55 


59 

60 

2900 

3000 

61 


62 

63 

64 


WRITE {6, 
WRITE (6 
WRITE (6, 
WRITE (6, 
WRITE (6, 
FORMAT (1 
FORMAT ( 1 
FORMAT (1 
FORMAT ( 1 
FORMAT (1 
FORMAT ( 1 


2400) (B(N),N 

,2500) (C(N), 

2600) (DLZ(N) 
2700) (XL (N) , 
2800) ( XM ( N ) , 

H , 1 HA/70 ( 1 H 
H , 1 HB/70 ( 1 H 
H , 1 HC/70 ( 1 H 
H , 1 HD/70 ( 1 H 
H , 1 HL/70 ( 1 H 
H ,1HM/70(1H 


= 1 , N XY) 

N= 1 , NX Y) 

, N= 1 , NXY) 

N= 1 , NX Y) 

N= 1 , N X Y) 

, 10D13. 5/) ) 
, 10D13.5/) ) 
.10D13.5/}) 
, 10D13.5/) ) 
, 10D13.5/) ) 
, 10D13.5/) ) 


CONTINUE 
PMIN=. IE-20 
PRMX=. ID-20 
DO 55 J= 1 ,N Y 
DO 54 1=1, NX 
N=I + NX* (J-1) 

JTWIN=N Y 1- J 

IF(J.EQ. NY) YG (I) =PR.(N) 

P2(I,J,1)=PR(N) 

DUMM (I,J,1)=PR(N) 

DUMM (I, JTWIN , 1) =PR (N) 

DUMM (I,J,2)=H (N) 

DUMM (I, JTWIN, 2) =H (N) 

IF (J. NE. J MIN .OR. PR (N) . LT. PM IN) GD TO 57 
PMIN = PR (N) 

NMIN=N 

CONTINUE 

IF (JTWIN. NE. UMAX .OR. PR ( N) . LT . PM IN) GO TO 58 
PMIN=PR (N) 

NU=N 

CONTINUE 

IF (PR (N) . LT. PRMX) GO TO 54 
PRMX=PR (N) 

PRMXS = PR (N) 

NHOLD=N 
XPK=XG (I) 

P2 (I, J, 2) =H(N) 

CONTINUE 
KMAX=NXY- 1 
HMAX=H (KM AX) 

IF (K (LMAX) .LT.HMAX) H MAX=H (LMAX) 

IF (H (LU) . LT. HMAX) HMAX=H (LU) 

IF (IT (K 1 ) .LT.HMAX) HM A X = H { K 1 ) 

IF (PMIN. LT.PR (KMAX) ) PMI N= PR (K M AX) 

IF (PM IN. LT. PR ( K 1 ) ) PMIN=PF(K1) 

CONTINUE 

CONTINUE 

WRITE (7,2900) ( ( (P2 (I, J, K) , 1= 1 , NX) , J= 1, » Y) , K - 1, 2) 

FORMAT (8 D 1 0. 5 ) 

WRITE (6, 3000) PRi X , NHOLD 

FORMAT ( 1 H ,6HPRMX = ,D16.5,10X,7HNHOLD=,I10) 
CONTINUE 

CALL FSETUP(NX,NYMl,XG,YGC,DUMM) 

CALL CTRSBR(NX,NYM1,NXY,NY,X3,YGC,DUMM,XMIN, YMIN) 
CONTINUE 

CALL GRAFNO(IN,X1AX,XG,YG,NX,INY) 

CONTINUE 
CALL TERM 
CONTINUE 
STOP 
END 
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TABLE I. - ESTIMATE OF PERCENT LOSS OF LOAD CAPACITY FOR 


VARIOUS GEOMETRIES AND INTEGRATION DOMAINS 


Minimum 

film 

thickness, 

H o 

Radius 

ratio, 

a 

Inlet parameter 

Load capacity 

Loss of load 
capacity, 
m, 

percent 

Abcissa 

-x E • 

Ordinate 

±y e 

Finite 

domain, 

w(N) 

Infinite 

domain, 

W oo(N) 

mm 

36.5361 

1. 

00 

7. 

00 

1.5162 

1. 5237 

0.49 

■EsSi 

1.0000 



1.00 

.1532 

.1540 

.51 

10" 4 

36,5361 



7.00 

.4744 

.4818 

1.55 


25.2897 





.3917 

.3977 

1.50 


15.8434 





.3055 

.3100 

1.46 


11.8266 





.2604 

.2642 

1.45 


8.3033 



4. 

00 

.2133 

.2165 

1.50 


5.2073 



4.00 

.1635 

. 1659 

1.46 


3.9983 



2. 

00 

. 1369 

.1391 

1.61 


2 . 842 8 



2. 

00 

.1091 

. 1109 

1.54 


2.3213 



2.00 

.0946 

.0961 

1.51 


1.3978 



1. 

34 

.0640 

.0650 

1.56 


1.0000 



1.00 

.0479 

.0487 

1.61 
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TABLE n. - DATA SHOWING EFFECT OF GEOMETRY ON 


MINIMUM FILM THICKNESS 


[Average surface velocity in x -direction, u, 10 cm /sec; fluid vis- 
cosity, Vq, 0. 411x10 ^ N* sec/cm 2 ; effective radius of curvature, 
R x , 1 . 11125 cm. ] 


Radius 

ratio, 

a 

Load 

capacity, 

w(N) 

Reduced 
hydro - 
dynamic 
lift, 

L 

Minimum film thickness 

Difference be- 
tween H. 

min 

““ H min , 

D, 

percent, 

Inputted 

value, 

H i»in 

Calculated 

from 

equation, 

n. 

min 

36.5361 

1.8450 

1.9020 

10 

-5 

a 0. 9785X10 ' 5 

-2. 14 


1. 8392 

1. 8960 

10 

-5 

.9789 

-2.11 


.5749 

1.8741 

10 

-4 

a l. 0079X10" 4 

.79 


.5704 

1. 8595 



1. 0075 

.75 

25.2897 

.4738 

1. 8714 



a l. 0075 

.75 

25.2897 

.4701 

1. 8568 



1. 0072 

.72 

15.8434 

.3681 

1. 8652 



a l. 0078 

.78 

15. 8434 

.3653 

1. 8510 



1. 0070 

.70 

11.8266 

.3128 

1. 8592 



1.0079 

.79 

11.8266 

.3104 

1. 8449 



1. 0072 

.72 

8. 3033 

.2571 

1.8658 



a . 9919 

-.81 

8. 3033 

.2552 

1. 8520 



.9907 

-.93 

5. 2973 

. 1955 

1. 8501 



a . 9900 

-1.00 

5. 2973 

. 1940 

1. 8360 



.9894 

-1.06 

3.9983 

. 1627 

1. 8367 



a . 9899 

-1.00 

3.9983 

. 1615 

1. 8232 



.9887 

-1.13 

2.8428 

. 1284 

1. 8206 



a . 9880 

-1.20 

2.8428 

. 1274 

1. 8065 



.9876 

-1.24 

2.3213 

. 1105 

1. 8074 



a . 9878 

-1.22 

2.3213 

. 1097 

1. 7943 



.9863 

-1. 37 

1.3978 

.0727 

1.7543 

1 


a l. 0061 

.61 

1. 3978 

.0720 

1.7384 



1. 0094 

.94 

1. 0000 

.0537 

1.7305 



a l. 0135 

1.35 


.0531 

1.7121 



1. 0200 

2.00 


. 1719 

1.7533 

10 

-5 

a . 9891 

-1.09 


. 1713 

1.7472 

10 

-5 

.9901 

-. 99 


Calculated by using the parabolic film assumption in theory. 




TABLE III. - EFFECT OF 


FILM THICKNESS ON 
REDUCED HYDRO- 


DYNAMIC LIFT 


Dimensionless 

Effect of H min 

minimum film 

on reduced 

thickness, 

hydrodynamic 

H . 

lift, 

mm 




10 A 

0.992 

10 ’ 5 

.997 



(a-2) x = 0 plane. 

(a) Two rigid solids separated by a lubricant film. 



(b-l)y = 0r (b-2) x *0 plane. 

(b) Equ ivalent system of a rigid solid near a plane separated by a lubricant film. 
Figure 1. - Contact geometry. 
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(c) Solution using Reynolds boundary conditions. 
Figure 2. -Effect of boundary conditions. 





Reduced hydrodynamic 




Figure 5. - Pressure profiles along rolling direction for two film thicknesses. 
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